Quantum Topological Excitations: from the Sawtooth Lattice to the Heisenberg Chain 
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The recently elucidated structure of the delafossite 
YCUO2.5 reveals a Cu-O network with nearly independent 
A chains having different interactions between the s — 1/2 
spins. Motivated by this result, we study the A chain for 
various ratios Jbb/Jbv of the base-base and base- vertex in- 
teractions. By exact diagonalization and extrapolation, we 
show that the elementary excitation spectrum, which (within 
numerical error) is the same for total spins Stot = and 1, 
has a gap only in the interval 0.4874(1) < Jbb/Jbv < 1.53(1). 
The gap is dispersionless for Jbb/ Jbv = 1, but has increasing 
fc-dependence as Jbb /Jbv moves away from unity, related to 
the instability of dimers in the ground state. 
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There is currently great interest in the topological exci- 
tations of quantum antiferromagnetic (AF) spin systems. 
The isotropic s = 1/2 Heisenberg chain has spinon exci- 
tations, while the anisotropic chain in its ordered phase 
has localized domain walls between opposite AF regions, 
both widely discussed in the literature. Later on, Shas- 
try and Sutherland [Q introduced another class of exci- 
tations: topological quantum solitons separating differ- 
ent regions of broken translational symmetry. The typ- 
ical example is the symmetric zigzag spin ladder, first 
addressed by Majumdar and Ghosh || (MG), in which 
nearest-neighbor (NN) triangles sharing a base site are 
also vertex-vertex coupled. Its lowest-energy excitations 
are so-called kinks (K) and antikinks (AK), both with 
similar characteristics and giving rise to a finite gap 
AE sa 0.234Ji when the interaction J2 between next- 
nearest- neighbors (NNN) is half that of NN spins Ji |§ . 

Then attention turned to the sawtooth or A chain, 
which consists of coupled s = 1/2 Heisenberg spins form- 
ing triangles aligned in a chain with a common base site, 
like the zigzag ladder but without the vertex- vertex cou- 
pling. The A chain has unique properties. In fact, nu- 
merical 0| and variational g studies of this lattice, all 
with bonds having the same interaction, have shown a re- 
markable feature: the K-AK symmetry of the MG model 
is broken here, yielding however a similar dispersionless 
small gap for the low-lying excitation modes. Very re- 
cently, in a study of an asymmetric zigzag ladder with 
NNN alternating interactions, that is, with J2 ± 6, Chen 
et al. H have discussed the crossover from the MG model 
(J2 = Ji/2, 8 = 0) to the A chain with equal interactions 
(6 = J 2 = Ji/2). 



Theoretical interest in the sawtooth lattice itself is now 
enhanced by new experimental results. In fact, the over- 
doped i?Cu02+i (R = Y, La, etc.) delafossite com- 
pounds, first synthesized by Cava et al. 0, opened up 
new possibilities for studying hexagonal Cu planes with 
AF interactions between the Cu 2+ ions. Depending on 
the O-doping, different s = 1/2 effective magnetic lat- 
tices are obtained, although with much weaker interac- 
tions than high-T c systems, which have comparable bond 
lengths but 180° Cu-O-Cu angles. Studies §,§ of the 
diluted kagome lattices of the x = 0.66 case predicted in- 
teresting properties. The detailed structure of YCUO2.5 
has just been elucidated [^0|JTl| l and, as suggested earlier 
Ijjjj, this compound appears as a physical realization of 
the sawtooth lattice (Fig. |l|). Recent synthesis has 
succeeded in obtaining an orthorhombic 2H single phase, 
in which the additional x = 0.5 oxygen ions are located 
at the center of alternating sets of triangles, providing 
super-exchange paths between s = 1/2 spins on nearly 
independent A chains. However, careful consideration 
of the measured angles and distances indicates different 
interactions between two spins on the base Jbb and be- 
tween the base- vertex NN spins Jb v of the triangles, most 
probably with Jbb < Jbv, owing to the smaller Cu-O-Cu 
angle ]lc| ]. While the case Jbb = Jbv has been quite ex- 
tensively studied theoretically |^|||, to the best of our 
knowledge the case Jbb 7^ Jbv has never been considered 
before. 

Therefore, we here analyze the sawtooth lattice for var- 
ious ratios Jbb/ Jbv of these AF couplings. The Hamilto- 
nian is given by 
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where N is the number of triangles (2N spins) in the 
chain, and Sj is the spin-1/2 operator at site i. Note that 
when Jbb = we retrieve the Heisenberg chain with 2A^ 
sites, while when Jb v = we have a Heisenberg chain 
with TV sites and N disconnected spins. It is therefore 
interesting to study the evolution of the elementary exci- 
tations from the sawtooth lattice with equal interactions 
Jbb = Jbv, hi which K-AK pairs yield a dispersionless 
small gap QJ^] of the same order as that in the symmet- 
ric zigzag ladder, to the isotropic Heisenberg chain, in 
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which pairs of spinons exhibit a strongly dispersive spec- 
trum JT2| with excitation energy varying between zero 
and irJ/2 (J is the coupling between adjacent spins). 

To test our numerical procedure, which includes larger 
clusters than before (up to 12 triangles), and to make 
contact with previous theory, we first briefly reconsider, 
and extend, existing numerical results for the A chain 
with Jbb = Jbv = J. In this case, Eq. (|l|) has two degen- 
erate ground states with N dimers pM. These ground 
states may be written as states in which each spin on the 
base of a triangle forms a singlet either with the follow- 
ing vertex spin (right, R-dimer state) or with the previous 
vertex spin (left, L-dimer state), that is, 

N 

|GS,R)=IJ[2<-1 ) 2<], (2) 

i=l 
N 

|GS,L) =JJ[2t,2t + l] 1 (3) 

i=l 

where = (|a</3j) - \0iaj))/V2, with a t ($) denoting 
the states with sf = 1/2 (—1/2) at the site i. These two 
states are linearly independent and become orthogonal 
for N — > oo. The existence of an excitation gap was also 
rigorously proved jH^]. The elementary excitations have 
been shown |||| to be well-separated K-AK-type domain 
walls separating regions of R-dimers and L-dimers. A K 
has a dimer in its triangle, while an AK does not. Cu- 
riously, they have very different characteristics in this 
system. A K has no excitation energy and is localized, 
but an AK propagates with kinetic energy within a region 
bounded by kinks. As a consequence of the former prop- 
erty the low-lying excitation spectrum is dispersionless, 
and owing to the second one, the gap value is consider- 
ably reduced compared to the energy of a trivial triplet 
state replacing a singlet dimer of the ground state. 

We diagonalize the spin Hamiltonian H in Eq. ([]]) by 
the Lanczos algorithm, using periodic boundary condi- 
tions with S2Ar+i identified with si . All sizes from N = 4- 
12 triangles are calculated. The values of the excitation 
energies obtained are then extrapolated to the N — > oo 
limit by assuming the finite-size error term to be a poly- 
nomial in 1/N, whose coefficients are determined by fit- 
ting. In Table | we give the gaps found for wavevectors 
k = 0, 7r/2, and tt for excitations with total spin Stot = 0, 
1, and 2, for the case Jbb = Jbv = J- (Note that the 
wavevector k of an eigenstate \tp) is here defined such 
that T n \ip) = e lkn \ip), where T n is the translation oper- 
ator by n triangles or 2n spins.) The quoted numerical 
error, which arises entirely from the extrapolation, is de- 
fined (here and later) to be twice the change in the result 
upon discarding the data of the largest system calculated 
(12 triangles) and repeating the extrapolation. In this 
case, the leading finite-size effect goes as 1/N 2 . Table ffl 
confirms the gap to be dispersionless, within numerical 
error. Our most precise value of the gap, for k = and 



Stot = lj is AE 1 = 0.2155(3)J, confirming the estimate 
0.215 by Nakamura and Kubo 0], and a little smaller 
than the value 0.219 expected from the variational cal- 
culation by Sen et al. 0] for an AK configuration spread 
over 5 sites. We find that, within numerical error, the 
Stot — low-lying excitations become degenerate with 
the Stot = 1 gap as A — + oo. Furthermore, the spectrum 
for Stot = 2 appears also to be dispersionless, but about 
twice that for Stot = 1 or 0, as shown in Table ffl. This 
new result is contrary to the speculation in Ref. H] that 
the excitation energies for higher spins might converge to 
the same limiting value as N — > oo. 

Figure 2 shows the evolution of the low-lying triplet 
(Stot = 1) excitation spectrum for Jbb/ Jbv < 1- As 
Jbb decreases, the triplet excitation energy decreases at 
k = until it vanishes near Jbb /Jbv ~ 0.5, while for 
k = 7r it goes up. Progressively a stronger fc-dispersion 
appears, yielding for Jbb = the famous lower bound- 
ary expression for the continuum of excited triplet states 
for the isotropic s = 1/2 Hciscnbcrg chain, A_E^(fc) = 
(7r/2) Jbv I sin/c/2| ]l2| (here rewritten keeping the defini- 
tion of k for our A chain) . 

For Jbb = Jbv the ground state is doubly degenerate 
for all N. But for Jbb < Jbv and N finite, these states 
are split into a ground state with spin-parity + and a 
first excited singlet state with spin-parity - . However, 
we find that the latter approaches the ground state with 
asymptotic dependence 1/N as as N — > oo, restoring the 
double degeneracy. 

An accurate determination of the critical ratio 
Jbb/ Jbv < 1 required to produce a triplet gap may be 
made by the method proposed in Ref. 01|. The idea 
is to map this problem to a continuum field theory and 
take into account the fact that in a fermion system with a 
fixed number of particles umklapp scattering is the only 
interaction that splits the degeneracy of the two lowest 
excited states. Therefore, the difference of their ener- 
gies provides a precise measure of the umklapp processes, 
which vanish at the critical interaction ratio. For the A 
chain, with N finite and Jbb < Jbv, the first two excited 
states have spin-parity _ and 1 + and undergo a level- 
crossing near Jbb/ Jbv ~ 0.5. As shown in Fig. |], the 
value of Jbb / Jbv at the level crossing follows a polyno- 
mial in 1/N 2 , and its extrapolation to N — > oo yields 
(Jbb/Jbv)crit = 0.4874(1). 

The ground state when Jbb = Jbv consists of spin sin- 
glets between NN spins, the R-dimers or L-dimers of Eqs. 
(^) and (^). To study the evolution of the ground state |0) 
for Jbb < Jbv j we plot in Fig. H the dimerization fraction, 
defined as £> frac - |(0|GS,R)p + |(0|GS,L)| 2 g|. D frac 
drops from its value of unity at Jbb = Jbv toward zero 
as Jbb/ Jbv decreases, showing that the gap reduction is 
related to increasing fluctuations of the dimer state. Fig. 
^ shows how important the extrapolation to ./V — > oo is in 
this calculation; we found that an Aitkens extrapolation 
gave very stable results. 
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Turning to the case Jbb/Jbv > 1, we show in Fig. || 
how the Stot = 1 gap dispersion curves are modified. 
The minimum gap is now found for k — n and decreases 
with increasing interaction ratio. On the other hand, 
for large interaction ratios the low-lying states become 
nearly degenerate. This is easy to understand: we are 
again approaching the Heisenberg chain, though just for 
the N sites on the bases of the triangles, while the re- 
maining N spins on the vertex are only loosely coupled, 
leading to a complex of 2 N nearly degenerate states. The 
evaluation of the critical ratio for the closure of the gap 
also becomes more complicated (details will be given else- 
where) , but an approach related to that used above gives 
( Jbb/Jbv)crit = 1-53(1). 

Figure 6 summarizes our main numerical results for 
the low-lying excitations of the sawtooth lattice. A fi- 
nite gap is found only for interaction ratios within the 
interval 0.4874(1) < J bb /Jbv < 1.53(1). Thus, the 
curve is not completely symmetric around Jbb/Jbv = 1- 
While our above discussion of the excitation gap has 
been explicitly just for Stot = 1 excitations, we note 
that we have also been able to calculate many features 
of the spectra for S tot = with a numerical accuracy 
of better than a few percent, finding agreement with 
the Stot = 1 values in all cases. Examples: the exci- 
tation energy at k = tt for < Jbb/Jbv ^£ 1.5 and at 
k = for 0.9 < Jbb/Jbv < 1.0. Also, for S tot = 
gap closure we find ( Jbb/ Jbv)crit = 1-51(3) compared to 
(Jbb/Jbv)crit = 1-53(1) for S to t = 1- Further details 
will be published elsewhere. This provides strong nu- 
merical evidence that the lowest excitation spectra are 
in fact four- fold degenerate (in the limit N — > oo) for all 
< Jbb/Jbv i5 1-5, thus generalizing the known result 
for the isotropic s = 1/2 Heisenberg chain (Jbb = 0) |L7|. 

To conclude, we have shown that different interactions 
along the base and vertex bonds in the sawtooth lattice 
yield a fc-dependent elementary excitation spectrum, the 
effective gap being progressively reduced with increasing 
instability of the dimer ground state. The recent crys- 
tallographic study of YCUO2.5 gives values for the three 
sides of the triangles and for the Cu-O-Cu angles. The 
consideration of these parameters and of the different co- 
ordination of the Cu ions (see Fig. @ and Ref. jl(j) justi- 
fies taking the same base-vertex interaction Jb v for both 
sides of the triangle, but a different Jbb along the base, 
as done here. The effect of a stronger AF interaction in 
just one of the vertex bonds is quite clear: it will lift the 
degeneracy of the gapped ground state, stabilizing the 
dimer singlet ordering corresponding to that direction. 
Preliminary susceptibility measurements on YCUO2.5 by 
Garlea and co-workers [jnj yield a low-temperature gap 
~ 17 K, which would imply AF interactions of ~80 K if 
one assumes Jbb = Jbv, a very reasonable value for this 
system. Now that good samples are available, measure- 
ments using more appropriate methods are in progress 
to distinguish the different interactions and to obtain a 



TABLE I. Values of the lowest-energy excitations of the 
sawtooth lattice for Jbb = Jbv with total spin Stot and 
wavevector k, after extrapolation to N — » 00. 



"Stot 


gap (k = 0) 


gap (k = tt/2) 


gap (k = 7r) 





0.2153(8) 


0.22(1) 


0.216(2) 


1 


0.2155(3) 


0.214(10) 


0.216(2) 


2 


0.46(1) 


0.49(8) 


0.46(6) 



precise value of the gap, if one exists (depending on the 
interaction ratio). It will be worthwhile to synthesize 
single crystals to allow the study of the dispersion of the 
excitation spectrum. 



* E-mail: sblundell@cea.fr (SAB), nunezreg@polycnrs- 
gre.fr (MDNR) 
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FIG. 1. Sawtooth chains in the triangular Cu planes of the 
delafossite YCUO2.5. The extra O ions (white) for x — 0.5 
are located at the center of particular triangles of Cu ions, 
creating AF super-exchange only within these triangles. This 
gives nearly independent A chains, indicated by thin black 
lines. While Cul (black) adopts tetrahedral coordination with 
two O ions in this plane (and two other O out of it), Cu2 (gray 
circles) adopts triangular coordination with just one O in this 
plane. The angles and distances Jl^] suggest a weaker Jbb 
interaction between Cul-O-Cul bonds (bases of the triangles) 
than for the Cul-0-Cu2 base- vertex bond Jb v - 



FIG. 2. Low-lying excitation spectra of the sawtooth lat- 
tice for Jbb /Jbv < 1. The dashed line corresponds to the 
Jbb = Jbv case, while the Jbb = dotted-dashed curve fol- 
lows the dispersion of the isotropic Heisenberg chain (using 
here the definition of k for the A chain). 



FIG. 3. Variation of the critical interaction ratio for gap 
closure when Jbb /Jbv < 1, as a function of the number of 
triangles N. The polynomial extrapolation in 1/N 2 yields 
the best estimation of this critical value. 



FIG. 4. Dimerization fraction (see text) vs. the interaction 
ratio. The calculations have been done for up to 12 triangles. 
The importance of the extrapolation is shown. 



FIG. 5. Dispersion curve for the gap to Stot = 1 excited 
states of the sawtooth lattice with Jbb/Jbv > 1, compared to 
the Jbb/Jbv = 1 case. 



FIG. 6. Gap to the lowest Stot = 1 excited states of the 
sawtooth lattice vs. Jbb/Jbv, after extrapolation to N — > 00. 
The error is greater for Jbb/Jbv > 1 because states with k — ir 
exist only for even N, giving less points in the extrapolation. 
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